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ABSTRACT 


A derivation  is  given  of  the  nonlinear,  small-disturbance  equations 
governing  three-dimensional  transonic  flow  through  a fan  or  compressor  rotor. 

1'hese  equations  represent  the  counterpart,  for  turbomachinery  flows,  of  the 
small-disturbance  equations  appropriate  to  an  isolated  airfoil.  Thus  they 
facilitate  the  application  to  turbomachinery  flows  of  the  many  numerical 
solution  methods  developed  in  recent  years  for  isolated-airfoil  problems. 

Boundary  conditions  for  design  and  off-design  conditions  are  formulated, 
at  a level  of  approximation  consistent  with  the  small-disturbance  field  equations. 
An  order-of-magnitude  analysis  is  presented,  which  reveals  the  transonic  simi- 
larity parameters  for  flow  in  blade  rows  of  low  and  high  solidity. 

The  validity  of  a rectilinear  transonic  shear  flow  as  a representation 
of  the  flow  through  a fan  or  rotor  is  also  discussed. 
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INTRODUCTION 


Axial-flow  compressors  have  been  operated  for  many  years  in  the 
transonic  range,  i.e.,  such  that  the  resultant  of  the  axial  approach  velocity 
of  the  flow  and  the  circumferential  velocity  of  the  blades  os  sonic  somewhere 
between  the  hub  and  the  tip.  Under  these  circumstances,  an  observer  riding 
on  one  of  the  blades  sees  a mixed  flow,  in  which  a subsonic  region  near  the 
hub  exists  side-by-side  with  a supersonic  region  near  the  tip.  These  two 
regions  interact,  making  three-dimensionality  an  essential  feature  of  the  flow. 

The  addition  of  this  feature  to  the  complexities  already  inherent  in  the  transonic 
character  of  the  flow  mak*.s  the  problem  a very  formidable  one. 

Despite  these  complications,  transonic  compressors  have  given  satis- 
factory performance,  and  in  recent  years  the  conflicting  demands  of  increased 
compression  per  stage,  with  decreased  weight  and  noise  generation,  have 
continued  to  place  great  emphasis  on  the  flow  field  problems  of  these  devices, 
which  are  inherently  three-dimensional  and  nonlinear.^ 

With  one  or  two  exceptions,  most  of  the  existing  studies  of  these  flow 
fields  have  separated  the  three-dimensionality  and  the  nonlinearity.  Thus, 
there  are  extensive  studies  in  which  the  nonlinearity  is  accounted  for,  but  in 
a two-dimensional  approximation  (see,  for  example,  References  2 and  3)  and 
a lesser  number  in  which  the  three-dimensionality  is  studied,  but  in  a linear 
approximation  (References  4 and  5 are  typical).  The  nonlinear,  two-dimensional 
solutions  are  often  applied  in  an  iterative  manner,  starting  with  an  axisymmetric 
calculation,  to  determine  a stream  surface  location  in  a meridional  plane.  Then 
the  flow  in  the  blade-to-blade  plane  is  calculated.  The  success  of  the  latter 
calculation  is  seriously  affected  by  the  presence  of  shock  waves  in  the  field; 
the  problem  of  locating  them,  and  of  adjusting  their  locations  on  subsequent 
iterations  has  been  found  to  be  very  difficult.  Thus,  the  three-dimensionality 
and  the  nonlinearity  become  especially  hard  to  handle,  even  in  an  iterative  way, 
when  transonic  flow  conditions  occur. 
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In  contrast,  many  calculations  of  the  transonic  flow  over  isolated  airfoils 
have  appeared  in  recent  years.7  **  These  methods  have  reached  a stage  of 
development  that  allows  three-dimensionality  and  nonlinearity  to  be  included 
simultaneously.  However,  it  appears  that  most  of  these  significant  advances 
have  not  yet  been  successfully  carried  over  to  the  calculation  of  turbomachinery 
flows . 

The  purpose  of  the  present  study  is  to  adapt  the  isolated-airfoil 
computational  techniques  to  the  problems  of  turbomachinery  flows.  As  the  first 
step  in  this  process,  the  appropriate  nonlinear  small-disturbance  equations  are 
derived  in  this  paper.  These  equations  are  presented  in  a form  that  facilitates 
application  of  virtually  all  of  the  isolated-airfoil  numerical  techniques.  The 
successful  adaptation  of  these  techniques  to  the  conditions  of  flow  in  a rotating 
blade  row  is  a logical  first  step  in  the  exploitation  of  the  whole  range  of 
methods  for  solving  the  fully  nonlinear  problem. 

The  first  section  below  retraces  the  derivation  of  the  basic  equations, 
paying  careful  attention  to  the  conditions  under  which  the  flow  seen  by  a 
blade-fixed  observer  may  be  considered  steady.  The  results  of  this  section 
shed  light  on  the  approximation  implicit  in  the  use  of  a rectilinear  transonic 
shear  flow  as  a model  of  the  flow  in  a transonic  compressor  blade  row. 

The  second  section  takes  up  the  special  case  in  which  the  flow  in 
blade-fixed  coordinates  consists  of  small  disturbances  from  a helical  stream. 

The  leading  nonlinear  terms  which  govern  the  transonic  character  of  the  flow 
are  retained,  in  a fully  three-dimensional  treatment. 

Since  the  problem  contains  several  small  parameters  (including  the 
departure  of  the  local  Mach  number  from  unity  and  also  the  geometric  and  loading 
parameters  that  produce  the  small  disturbances  from  helical  flow)  an  order-of- 
magnitude  analysis  is  needed,  to  identify  the  relative  sizes  of  the  parameters 
that  may  consistently  be  retained.  This  analysis  is  presented  in  Section  III. 

The  concluding  section  reviews  some  of  the  implications  which  the  present  results 
will  have  for  later  numerical  work,  and  extensions  to  more  fully  nonlinear 
approximations . 
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I.  THE  BASIC  EQUATIONS 

The  purpose  of  the  present  work  is  to  facilitate  the  application,  to 
turbomachinery  flowfields,  of  a number  of  numerical  methods  which  were  developed 
originally  for  external  flows  over  isolated  bodies.  Thus,  although  the  equations 
of  motion  for  turbomachinery  flows  are  well  known,  nevertheless  there  is  a need 
to  retrace  their  derivation,  in  order  to  display  the  forms  analogous  to  those  that 
have  been  used  in  numerical  studies  of  external  flows. 

Coordinate  Systems 

Several  coordinate  systems  can  be  used  to  describe  the  flow  in  a 
compressor  blade  row.  The  first,  called  "space-fixed”,  is  one  in  which  the  rotor 
advances  and  rotates.  The  flow  in  this  coordinate  system  is  always  time-dependent. 
The  subscript  ( ) is  used  in  this  paper  to  designate  these  coordinates.  (See 

Figure  1.) 

The  second  system  can  be  thought  of  as  being  fixed  in  the  engine  mount. 

In  these  coordinates,  the  blade  row  rotates,  but  does  not  advance.  Following 
Wu*5,  these  coordinates  are  referred  to  as  "absolute".  They  are  designated  here 
by  the  subscript  ( ) . The  flow  is  unsteady  in  these  coordinates.  (See  Figure  2.) 

A third  system  of  coordinates,  introduced  below,  has  an  origin  fixed  on 
one  of  the  blades.  This  system,  referred  to  by  Wu  as  "relative",  is  denoted  in  this 
paper  by  the  term  "blade-fixed",  and  by  coordinate  symbols  with  no  subscript. 

Under  certain  conditions,  the  flow  in  these  coordinates  may  be  regarded  as  steady. 

Equations  of  Motion 

In  the  absolute  coordinates,  the  component  forms  of  the  continuity  and 
momentum  equations  are: 
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The  quantities  Vj_  , , 6^  arc  cylindrical  coordinates,  with  V 2.  along  the 

axis  of  the  blade  row,  in  the  radial  direction,  and  measured  positive 

in  the  right-handed  sense  (see  Figure  1). 


Next,  these  equations  are  written  in  blade-fixed  coordinates,  using 
velocity  components  relative  to  the  blades.  The  relative  velocity  (V  is 
given  by: 

\J  = + tJg  * f 


where  uJg  is  the  angular  velocity  of  the  blade  row: 

A 

iJg  ~ 

A 

where  is  a unit  vector  in  the  y -direction.  Then 
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The  blade-fixed  coordinates  (see  Figure  3)  are  given  by 
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In  terms  of  these  variables,  the  continuity  and  momentum  equations  become: 
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Note  that  these  equations  are  unchanged  if  the  signs  of  ail  the  angular  variables 
are  reversed.  This  is  equivalent  (see  Figure  3)  to  reversing  the  sense  of  the 
angular  variable  and  the  angular  velocity  component,  and  to  interpreting  u)  as 
the  angular  velocity  of  the  flow  relative  to  the  blades: 


\fj9  S ^ ) & » - + Const  aJ  = - ujg 

If  there  are  no  flowfield  disturbances  which  vary  with  <9*,  in  the  absolute 
coordinates,  then  the  flow  in  the  blade-fixed  coordinates  may  be  taken  as  steady. 
This  means  that  there  can  be  no  rows  of  stators  or  inlet  guide  vanes  nearby,  nor 
any  non-axisymmetric  inlet  distortion.  Thus,  strictly  speaking,  the  rotor  must 
be  completely  isolated.  This  restriction  is  very  carefully  stated  by  IVu  (Refer- 
ence 15,  p.  12)  and  by  Vavra  (Reference  16,  pp.  111-112). 


In  actual  practice,  it  may  be  expected  that  the  steady-flow  approxi- 
mation in  blade-fixed  coordinates  will  have  acceptable  accuracy  whenever  the 
neighboring  stators  are  several  chord  lengths  away.  The  fact  that  this  approxi- 
mation has  been  used  routinely  for  many  years  (see,  for  example.  Reference  17) 
for  design  purposes,  even  at  very  small  rotor-stator  spacings,  is  evidence  that  the 
steady-flow  approximation  is  satisfactory  for  at  least  some  purposes.  In  the 
remainder  of  the  present  work,  this  approximation  will  be  adopted,  in  common  with 
the  bulk  of  existing  work  on  turbomachinery  flowfields.  Thus,  the  basic 
equations  used  are  Eqs . (8) -(11) , with  time  derivatives  set  equal  to  zero, 
replaced  by  & , and  —oJg  replaced  by  u>  . 


The  law  of  conservation  of  energy,  expressed  in  the  absolute  coordinates, 
states  that  the  change  in  total  enthalpy  per  unit  mass  /■)'"$,+  on  passing 
through  the  blade  row  is  equal  to  the  mechanical  work  done  by  the  blade  row  on 
the  fluid.  In  blade- fixed  coordinates , it  is  much  more  convenient  to  use  a 
quantity  X,  defined  as 

x = + 

15 


W ' 


(ojr) 


(13) 


i Z (12) 

Wui,J  has  shown  by  straightforward  manipulations  that  X satisfies  the  equation 

1)1  - , i ii 

Dt  Dt  % 31 

where  S is  the  entropy  per  unit  mass,  and  denotes  the  convective 

derivative.  Thus,  in  a flow  which  is  steady,  and  where  the  entropy  is  constant 
along  streamlines,  the  quantity  I will  also  be  constant  along  streamlines. 
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These  approximations  are  adopted  here.  Thus  using 
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The  conservation  of  quantity  X along  streamlines  becomes; 
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The  terms  on  the  right-hand  side  of  this  equation  are  considered  to  be 
evaluated  at  some  arbitrary  point  ( * , r , & ) in  the  flowfield,  while  those  on 
the  left-hand  side  are  evaluated  at  a related  point,  on  the  same  streamline,  far 
upstream  of  the  rotor.  Since  the  flow  at  infinity  is  assumed  not  to  vary  with 
x or  & , the  terms  on  the  left-hand  side  vary  at  most  with  f"co  . From 
Equation  (.6) , it  can  be  seen  that 

X ^40  = ■£  j (X<u  + O + ( U)  )’U  j 

Thus  Equation  (14)  can  be  written  as 
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In  the  present  work,  both  (2®  and  will  be  taken  as  constants.  It  should  be 

stressed,  however,  that  any  nonuniformities  present  in  the  flow  far  upstream  of 
the  blade  row  are  convected  along  streamlines.  For  example,  if  the  gas  tempera- 
ture far  upstream  had  a slight  radial  variation,  then  the  local  sound  speed  would 
vary  as  » 


x *•  d^oo  (fp*)  , v 
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where,  to  first  order,  the  streamline  displacement  is 
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where  the  integration  is  done  along  an  undisturbed  streamline. 
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A single  equation,  containing  only  velocity  components,  can  be 
derived  by  combining  the  continuity,  momentum,  and  energy  equations.  First, 
us  ing 

^ - cf-  ^ ^ A. V 


the  continuity  equation  is  recast  as 
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7 J?  ^7 


If  the  and  & components  of  the  momentum  equation  are  multiplied  by 

Wy  , 4/r  > and  VJfi  , respectively,  and  added,  the  resulting  equation  can  be 


used  to  eliminate  the  pressure  terms  in  the  equation  above,  in  favor  of 
velocity  components.  Similarly,  the  energy  equation  can  be  used  to  replace 
derivatives  of  C2  . The  result,  after  some  rearrangement,  is: 
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II.  SMALL-DISTURBANCE  EQUATIONS 


Far  upstream  of  the  rotor,  the  components  of  the  velocity  vector  W 
approach  the  axial  velocity  and  the  wheel  velocity  which  are 

the  components  of  the  helical  approach  flow  seen  by  a blade-fixed  observer. 
In  what  follows,  attention  is  restricted  to  blade  designs  which  are  thin, 
and  essentially  aligned  with  these  helical  paths,  so  that  only  small  distur- 
bances from  free-stream  conditions  will  occur.  These  small  disturbance 
quantities  are  defined  as: 


W*  = iXco4  > W V * v'  , ^<9  “ 


cjr  + c^r 


(18) 


These  expressions  are  now  substituted  into  Equation  (17),  and  the  resulting 
expression  is  expanded  into  terms  which  are  linear,  quadratic,  and  cubic  in 

X. 

the  small  quantities.  In  addition,  the  quantity  A is  written  as 

a u = #2  + ( Aa>) 


where  Equation  (15)  gives 
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After  adding  and  subtracting  certain  terms,  so  as  to  display  the  vorticity 
components  explicitly,  the  result  is: 
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In  deriving  this  equation,  it  has  been  assumed  that  the  entropy  is  constant 
along  streamlines.  Thus  the  result  is  not  completely  general,  although  it 
does  contain  all  the  possible  velocity  terms  - no  small  quantities  have  been 
neglected  yet. 

In  the  remainder  of  this  paper,  the  disturbance  quantities  will  be 
assumed  small,  and  the  leading  terms  that  must  be  kept  in  the  transonic 
range  will  be  sought.  The  assumptions  of  small  disturbances,  and  particle  - 
isentropic  flow  imply  that  the  maximum  relative  speed  of  the  blades  must  be 
only  slightly  greater  than  sonic  (maximum  relative  Mach  numbers  on  the 
order  of  1.2  to  1.3,  say),  so  that  any  shockwaves  generated  will  produce 
negligible  changes  in  entropy.  Consistent  with  this  assumption,  the  flow  is 
taken  to  be  irrotational.  Thus,  any  vorticity  produced  by  shear  forces  at  solid 
surfaces  is  assumed  to  be  confined  to  boundary  layers  and  wakes  of  negligible 
thickness.  The  net  effect  of  these  assumptions  is  to  restrict  the  applicability 
of  the  results  to  thin,  lightly  loaded  blades,  operating  at  tip  speeds  which 
are  slightly  supersonic.  In  addition,  this  analysis  shares,  with  the  bulk  of 
turbomachinery  literature,  the  steady-flow  assumption,  which  rules  out 
cases  where  the  blade  row  is  located  close  to  any  stator  rows. 

These  assumptions  impose  an  upper  limit  on  the  range  of  blade 
geometries  and  operating  conditions  where  the  present  approximation  will 
give  good  agreement  with  experiment.  Presumably,  any  design  which  is  to 
be  efficient  at  transonic  speeds  must  employ  thin,  lightly  loaded  blades.  In 
the  case  of  isolated  airfoils,  there  is  a large  range  of  thickness,  loading, 
and  Mach  number  where  small-disturbance  theory  compares  favorably  with 
experimental  data.  For  the  case  of  turbomachinery  flows,  the  extent  of  this 
range  remains  to  be  determined,  by  comparison  of  its  predictions  with 
appropriate  experiments. 

The  approximations  inherent  in  the  small-disturbance  approximation  are 

1 8 

further  clarified  by  the  work  of  Kerrebrock  , who  subdivided  the  disturbances 
into  pressure,  shear,  and  entropy  modes.  He  states  that  these  three  modes  are 
in  general  coupled  to  one  another  in  turbomachinery  flows,  in  contrast  to  their 
independent,  additive  behavior  in  external  flows.  The  present  approximation 
admits  only  the  pressure  mode.  Shear  and  entropy  disturbances,  as  noted  above, 
are  not  present  in  the  inlet  flow,  and  are  not  created  during  passage  through 
the  blade  row. 
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The  absence  of  entropy  disturbances  is  assured  (see  Equation  15)  by 
adopting  a uniform  inlet  temperature,  while  shear  disturbances  are  excluded 
(see  above)  by  allowing  only  constant  values  of  dZ^and  tJ  . It  is  important 
to  note  that  the  radial  variation  of  the  resultant  velocity  does  not  constitute 
a shear  flow.  This  point  is  discussed  further  below. 

The  leading  terms  in  Equations  (20),  which  dominate  in  the  transonic 
range,  must  now  be  found.  It  is  clear  that  the  terms  on  the  last  three  lines 
can  be  dropped,  as  well  as  the  quadratic  terms  in  the  first  bracket  on  the  third 
line.  In  order  to  determine  which  of  the  remaining  terms  should  be  kept,  it  is 
convenient  to  use  a coordinate  system  S > ^ H aligned  with  the  helical 
free- stream: 

V 
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It  should  be  noted  that  the  S , A coordinates  do  not  coincide  with  the  local 
streamline  direction,  but  only  with  the  orientation  of  the  free-stream  direction 
In  general,  the  local  streamlines  have  a first-order  departure  from  the  free- 

stream  direction. 


The  chain  rule  expressions  for  the  partial  derivatives  become: 
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Thus  Equation  (20)  can  be  written  as 
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Note  that  the  radial  derivative  has  been  left  untransformed.  These  equations 
are  now  written  in  terms  of  velocity  components  that  lie  in  a surface  of  con- 
stant radius,  and  are  directed  along  and  normal  to  the  free-stream  direction: 
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It  can  be  shown  that 


as  *w  r 
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Consequently,  the  basic  differential  equation  becomes: 


(27) 
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This  equation  can  be  simplified  by  neglecting  the  products  of  disturbance 
quantities.  However,  in  the  transonic  range,  the  quantity  V/V0  ~ <^tx>  is  itself 
a small  quantity,  which  multiplies  the  term  3s  . Thus,  the  nonlinear 

terms  which  contribute  to  this  coefficient  are  retained,  giving 


r y V.  . I ^ ^ f J.  — 

[\A/0  - Qo  + ^ 1 “**»  Lz*  & r J (29) 


It  is  significant  that,  in  the  small-disturbance  approximation,  the  radial 
velocity  component  does  not  contribute  to  the  transonic  nonlinearity. 

In  addition  to  W0  “ ®~ao  , other  small  parameters  related  to  the  blade 
thickness  and  loading  will  play  a role  in  the  problem,  and  an  order-of-magnitude 
analysis  is  required,  to  establish  the  relative  magnitudes  that  may  be  retained  in 
a consistent  analysis.  This  subject  is  taken  up  in  Section  III. 

Transonic  Shear-Flow  Approximation 

The  similarity  of  Equation  (29)  to  the  conventional  two-  and  three- 
dimensional  rectilinear  transonic  small-disturbance  equations  is  obvious.  However, 
there  are  also  some  very  marked  differences,  which  are  worthy  of  further  examina- 
tion, since  rectilinear  transonic  shear  flows  are  sometimes  suggested  as  being 
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representative  of  the  flow  in  a transonic  compressor  blade  row. 
rectilinear  equations  used  in  these  papers  do  not  contain  any  of  the  axisymmetric 
terms  in  the  divergence  of  the  velocity  vector;  thus  in  the  small -disturbance 
limit,  their  use  is  tantamount  to  rewriting  Equation  (29)  as 

+ §£)  ] - « 

2 A a,<-  h , * ) (3  0) 


( \tjl~  + <YVi) 


i?Hi  "*• 
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Here  5 , ^ , and  V)  are  regarded  as  rectilinear  coordinates,  and  the  approach 
velocity  W0  is  taken  to  be  a linear  function  of  tZ  , rather  than  the  square- 
root  function  corresponding  to  the  resultant  of  (J-to  and  o-lf.  The  correct 
form  of  Equation  (29)  is 

r * » , /_*,.»  . 1 f av-  ^c\  av- a m\  v-l 

57  '““(w  *55  )r  + 55  77  » 5?)  * rr° 

(31) 

Thus,  the  rectilinear  model  can  be  thought  of  as  neglecting  three  terms:  one 

of  them  is  the  term,  characteristic  of  the  divergence  in  an  axisymmetric 

flow,  while  the  other  two  are  related  directly  to  the  radial  variations  of  the 
distances  measured  along  and  normal  to  the  streamlines. 


Moreover,  the  use  of  a rectilinear  shear  flow  also  introduces  a shear 

disturbance  that  is  not  present  in  the  uniform-inlet  turbomachinery  case.  This 

21  22 

point  is  illustrated  in  the  papers  of  Namba  and  Inger  , which  start  from  the 
full  rectilinear  three-dimensional  equations,  and  then  proceed  to  the  equations 
for  small  disturbances  of  a transonic  shear  flow  profile.  Their  equations,  for 
adiabatic  flow,  show  an  additional  term  on  the  right-hand  side  of  Equation 
(30),  proportional  to  U * ^ °/<fx  > in  addition  to  the  other  discrepancies  noted 
above.  This  additional  term  does  not  arise  in  the  present  analysis,  unless 
the  approach  flow  C£<n  itself  is  taken  to  vary  with  V . The  point  to  be 
stressed  is  that  the  radial  velocity  variation  produced  as  the  resultant  of 
constant  vectors  and  is  not  correctly  modeled  by  a rectilinear 

shear  flow,  since  small  disturbances  of  the  latter  flow  lead  to  an  equation 
which  differs  from  the  one  appropriate  to  steady  flow  in  blade-fixed  co- 
ordinates, 

Velocity-Component  Equations 

Equation  (29)  can  now  be  transformed  back  to  the  cylindrical  co- 
ordinates, using: 
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The  result  is: 

1 1 (< V1)  ■ 6r+0Kto  (*  +/>^r) 


3?  / PZ 


^ 2<5> 


cJ-V- 


; 


/> 


& '-■ 
<*$> 


^6* 

^5 


7 


1 dt 


as 


+ ( 


'vVf7)=o 


(33) 


To  this  must  be  added  the  statement  that  any  two  of  the  three  vorticity  com- 
ponents are  zero: 
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Potential  Equation 


The  differential  equation  to  be  satisfied  by  the  velocity  potential  can 
now  be  found  directly  from  Equation  (33).  Let  <P  denote  the  potential  whose 
gradient  gives  the  dimensional  perturbation  velocity  components;  then  let 
dimensionless  perturbation  quantities  be  defined  as: 
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Then,  denoting  partial  derivative,  by  subscripts: 
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The  connection  between  this  equation  and  its  linearized  counterpart 
(References  4,  5,  23-29  ) can  be  seen  by  rewriting  it  as  follows,  where  the 
nonlinear  term  has  been  taken  to  the  right-hand  side: 

- (£*  + &»)  ( &Z2-  4 ”*■  ^&a)  (37) 

' / 
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Okurounmu  and  McCune  have  derived  an  equation  similar  to  (37), 
in  which  the  right-hand  side  is  replaced  by  A 4-S^,  where  A is  a 
constant,  whose  magnitude  can  be  estimated  from  the  rotor  geometry  and 

j*5 

operating  conditions.  Handling  of  the  nonlinear  term  by  the  introduction  of 

this  constant  has  been  a very  useful  approximation  in  the  treatment  of  other 

30 

transonic -flow  problems. 

A complete  potential  equation,  in  which  small  disturbances  are  not 
assumed,  and  is  not  approximated,  was  presented  by  Wu  (see  Reference  15, 
Equation  (26)).  It  can  be  checked  that  the  pre  sent  re  suits  are  consistent  with  his. 

Coordinate  Transformation 

The  cylindrical  coordinates  2^  P & have  two  features  that  limit  their 

usefulness  for  finite -difference  evaluations  of  Equations  (33),  (34)  and  (36). 

The  first  is  that  the  blade  surfaces,  where  boundary  conditions  are  to  be 

applied,  do  not  lie  on  constant  values  of  ? , , or  ^ . The  second  is  that 

the  axis  of  the  forward  Mach  cone  through  any  point  does  not  lie  parallel  to 
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any  of  these  axes.  The  results  of  South  and  Jameson  show  that  unwanted 
numerical  problems  may  occur  when  the  Mach-cone  axes  are  not  nearly  aligned 
with  a coordinate  direction. 

The  S , /?  , kl  coordinates  (Equation  (21))  would  take  care  of  the 
latter  problem;  however,  they  are  not  free  from  the  first  limitation, 
because  the  variation  of  S with  f'  causes  the  blade -surface  locations  to  be  a 

’"Their  result  (see  Eq.  14,  p.  1374  of  Ref.  26)  does  not  contain  the  quantity 
V on  the  right-hand  side  of  Eq.  37,  apparently  because  they  neglected 
the  contribution  from  the  term  involving  in  Eq.  20. 
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complicated  function  of  the  radius.  Moreover,  the  fact  that  the  quantities 

q in  the  chain-rule  derivatives  are  nonzero  leads  to 
a maior  increase  in  the  number  of  terms  in  the  differential  equations, 
especially  for  the  velocity  potential. 

An  alternate  set  of  coordinates  that  appears  to  circumvent  the  limita 
tions  cited  above  is  the  skew  set  using  the  % -coordinate  introduced  by 
Goldstein"^"  : 
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In  these  terms,  the  velocity-component  equations  become: 
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and  the  velocity  potential  satisfies: 
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The  -derivatives  in  this  equation  are  taken  at  constant  T and  j*  ; thus, 

they  are  aligned  with  the  relative  freestream  velocity  vector.  The  supersonic  - 

33,  34 

propeller  studies  of  Ordway  and  Hale  ’ have  shown  that  the  forward 
Mach  cones  through  any  point  are  centered  on  the  helical  streamlines,  and 
are  closely  approximated  by  a right  circular  cone,  in  the  vicinity  of  the 
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Boundary  Conditions  at  the  Blade  Surfaces 

The  Z , £ coordinates  have  the  advantage  of  convenience  in 

expressing  the  boundary  conditions  at  the  blade  surfaces.  If  there  are 
8 blades,  whose  chord  varies  with  \T  in  such  a way  as  to  have  a constant 
axial  projection  , then  they  are  located  at 


r = 

3 


■< s-'  ; 2 ' zr 


The  conditions  prescribed  at  the  blades  will  involve  either  the  blade  shape 
or  the  local  pressure;  these  are  related  to  the  velocities  by  (see  Reference  5, 
P.  9): 


)^/o_  \ 

'J-GD  / wjj 


(Here  ^ denotes  the  density  far  upstream).  The  quantities  Ur  and 
are  related  to  the  other  velocity  components  and  the  potential  by: 

= (~f)  (£--  *u)  = i § \ _ J—  £ \ 
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These  perturbation  quantities  will  appear  in  the  boundary  conditions 
in  various  combinations,  depending  on  how  the  problem  is  posed.  There  are 
at  least  three  cases  of  interest: 
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1.  ) The  first  case  might  be  referred  to  as  the  "off-design"  condition.  In 
it,  the  entire  blade  shape  is  given  in  terms  of  the  camberline  and  thickness 
distributions,  and  £ ( respectively: 

fY\,.  (*)  = ~ t (s)/a. 


- ('*■«+  x^/%  ; -t  ( ss  = -n«.  - %. 


Z,)  In  the  second  case,  which  might  be  called  the  "design"  condition,  the 
blade  thickness  distribution  f fi)  and  the  loading  distribution  on  the  blade 
■/j?LCs^  - ^(j'jare  given.  Thus,  the  boundary  conditions  will  involve  6/-^  , and 
differences  of  U.s  across  the  blade  surface. 

3.  ) A third  case  of  interest  is  that  in  which  the  pressure  distributions 
on  upper  and  lower  surfaces,  'j’JS)  and  "j^(s),  are  specified  (as  might  be 
dictated  by  boundary-layer  growth  considerations,  for  example).  The  problem 
is  then  to  determine  the  corresponding  blade  shape.  This  case  amounts  to 
a specification  of  Cls  on  each  blade  surface. 

In  all  of  the  above  cases,  the  boundary  conditions  are  applied  on  the 
line  ^ x rather  than  at  the  blade  surfaces.  This  linearization  of  the 

boundary  conditions  is  co.  sistent  with  the  order  retained  in  the  differential 
equations. 

At  the  hub  ( ) and  at  the  casing  (P~/\  )>  the  radial  velocity  is 

proportional  to  the  streamwise  slope  of  the  hub  and  casing  surfaces.  In  the 
present  work,  these  surfaces  are  assumed  to  be  cylindrical;  thus,  the  radial 
velocities  at  these  walls  are  zero: 


£ </>„')  ’ 5.  C/r)  =0 
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Periodicity  and  Farfield  Conditions 

The  flowfield  must  be  periodic  in  the  £»  -direction  at  fixed  -£  , 
with  period  "£,r  /£,  . Thus  the  perturbation  quantities,  such  as  velocity 

components  and  pressure,  must  display  this  periodicity.  However,  it  is  not 
necessary  that  the  velocity  potential  be  periodic. 

In  an  analytic  solution  that  is  constructed  by  an  eigenfunction 
expansion,  it  is  easy  to  satisfy  the  periodicity  requirement  by  proper  choice 
of  the  eigenfunctions.  For  the  present  considerations,  which  are  directed  at 
numerical  solutions,  this  procedure  is  not  available.  However,  an  alternate 
procedure  that  appears  to  be  suited  to  the  mixed  elliptic/hyperbolic  nature  of 
the  problem  is  to  require  that  the  perturbation  velocities  and  be 

periodic  on  the  helical  surfaces  in  which  the  blades  lie,  i.e., 

\ A ^ (*,/>,  o') 

(47) 

A 2F/<3)  = 


It  is  not  required  that  the  radial  perturbation  velocity  or  the  potential  be 
periodic  on  these  surfaces.  In  fact,  vortex  sheets  will  in  general  trail 
downstream  of  the  blade  surfaces,  and  both  the  potential  and  radial  velocity 
will  be  discontinuous  across  these  surfaces. 


The  perturbation  potential  is  defined  to  be  zero  far  upstream;  thus 
at  any  point  along  the  helical  surfaces  that  contain  the  blades,  the  potential 
is  given  by  2. 

* X ^ A°+) 

1 (48) 

r;i  * X 4>e  (e,A 

where  = 2^/5  Since  the  integrands  are  periodic  upstream  of  the  blades, 

it  follows  that  itself  is  also  periodic  there: 

4>  (49) 

If  the  integration  is  now  extended  onto  the  blades,  the  results  will  in  general 
be  different,  until  the  trailing  edge  is  reached;  the  potential  difference 
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existing  there  will  remain  unchanged  downstream: 


I 


4>(*, /»,<>*)  -4  (*,/>,  ~ 

^ 190  (50) 

This  jump  in  potential  is  proportional  to  the  normal  force  per  unit  span 

i by  integrating  the  pressure  difference  (pressure  side  minus  suction 
side)  along  the  blade  at  constant  radius,  it  is  found  that 

77  / 7 u*  (si) 

The  jump  in  potential  is  also  related  to  the  mean  turning  angle  downstream  of 
the  blades 


f <i>  d t ^ 


Zt 


In  addition,  the  jump  in  potential  is  related  to  the  total  pressure  ratio 

i 7 

across  the  blade  row.  In  the  full  nonlinear  treatment  , this  quantity 
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is  found  from  the  stagnation  temperature  ratio 

__  / Toi  \ 

\o\  “ ( To.  / 

which  in  turn  is  given  by  the  relation: 

Cp  (t01  - T’oi  ^ - «o  ( ~ r,  V©v ^ 

Here  the  subscripts  / and  £ refer  to  stations  upstream  and  downstream  of  the 
blades.  In  the  present  notation, 

Me  - W9-^r  =‘*r 

If  the  mean  value  calculated  above  is  used  for  US'  , if  7o/  is  expressed  in 
terms  of  the  inlet  Mach  number,  and  if  streamline  deflections  r*  - 1 7 are 
neglected,  then  the  leading  term  in  the  expression  for  the  total  pressure  ratio 


V = 


=.  \ + 


t. 

^ Moo 

„ Xzl  M 1 

~ f loo 


— 

Zt 


The  vortex  sheets  that  trail  downstream  of  the  blades  are  assumed 
to  lie  in  the  helical  surfaces.  The  jump  in  potential  across  this  surface 
must  be  taken  into  account  when  making  finite-difference  approximations  to 
derivatives  near  the  sheet.  One  of  the  manifestations  of  this  problem  involves 

the  evaluation  of  . Following  the  procedure  outlined  by  Ballhaus  and 

12  * * 

Bailey  , the  potential  is  divided  into  odd  and  even  parts: 


^ 4>0  4-4/ 


where 

= k [ + t-t.) 

*1  VM-4  (V, 

Along  ^ « V , the  quantity  4^  is  zero,  because  is  independent  of  ; 

thus  by  differentiation,  +».  is  also  zero  there.  Across  . Mv i and 

are  continuous,  hence  and  are  continuous  also.  Thus  <£•  is  zero  on 

, which  says  that  the  entire  variation  of  ^ is  in  the  even  part. 


■ 


Consequently,  along  , the  even  part  of  <^>  satisfies  the  original 

partial  differential  equation,  and  can  be  differenced  in  a straightforward  manner, 
while  the  odd  part  satisfies 

A°  , <+,°  4.  JL  cb°  * o 

But  it  is  also  true  that  Thus  the  odd  part  of 

is  found  by  a finite-difference  approximation  of 

l*' ) * -Jj?)  { (MV  * i (i4>V  \ <54) 

Farfield  Conditions 

At  large  distances  downstream  of  the  blade  row,  the  perturbation 
velocities  must  not  change  along  the  helical  directions,  i.e.,  there  must  be 
no  variation  with  2 at  constant  f*  and  ^ : 

A f * * ^ t 3 ^ 

Thus  the  perturbation  potential  must  be  of  the  form 


The  value  of  the  second  term  can  be  found  by  the  requirement  of  mass  conservation: 

>n  = ^ tx*,  yr  ( rT'-  r*)  = J 

r ^ Q 

If  the  density  and  axial  velocity  are  expanded  to  first  order  in  the  pertur- 
bation velocities,  the  expression  for  the  corresponding  mass-flow  perturbation, 
which  must  be  zero  is: 


J^r  J ^ ) U - />  i^oo  ur  j dr,  = O 

The  integral  involving  t/r  is  related  to  44  (see  Eq.  52),  and  the  quantity 
£ is  equal  to  the  constant  C , which  can  be  solved  for  as 

— ft 

c = r — / (5£ 

TT  ( /-M^)  (/r  ‘A1)  7h 

This  constant  is  proportional  to  the  static  pressure  rise  through  the  blade 
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row,  as  noted  by  Okorounum  and  McCune. 
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Finally,  a Kutta  condition  must  be  imposed  at  the  trailing  edges  of 
the  blades.  In  a full  nonlinear  treatment,  this  would  require  that  a stagnation 
point  be  located  at  a rounded  trailing  edge  or  at  a trailing  edge  having  a non- 
zero wedge  angle,  and,  at  a cusped  trailing  edge,  would  require  that  the  pressure 
be  continuous  across  the  trailing  edge.  In  the  small-disturbance  approximation, 
it  is  not  possible  to  represent  a stagnation  point,  and  thus,  the  solution 
becomes  singular  at  the  trailing  edge.  The  finite-difference  evaluation 
of  the  solution  smoothes  out  this  singular  behavior,  replacing  it  with  finite 
values  of  the  velocity  components.  In  the  present  work,  the  circulation 
is  adjusted  until  the  pressures  corresponding  to  these  finite  velocity  com- 
ponents are  continuous  across  the  trailing  edge. 
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III.  SIMILARITY  PARAMETERS 

The  equations  presented  in  the  previous  section  contain  the  leading 
nonlinear  terms  that  must  be  retained  when  the  difference  of  the  relative  Mach 
number  from  one  is  small.  The  specific  combination  of  terms  retained  implies 
an  assumption  that 

3-  = o(i-«Z) 

Wo 

However,  the  boundary  conditions  imposed  at  the  blade  surfaces  also  require 
that  U$  and  U , or  some  combination  of  them,  will  have  orders  of  magnitude 
that  are  dictated  by  the  small  parameters  that  measure  the  blade  geometry  and 
loading.  Because  these  components  are  connected  by  the  equations  of  motion, 
the  small  parameters  must  then  be  related  to  each  other,  in  order  that  the 
equations  be  consistently  ordered  in  the  limit  where  these  parameters  approach 
zero.  This  section  contains  a discussion  of  the  geometric  and  loading  parameters, 
and  their  relation  to  in  the  small -disturbance  limit. 

The  blade  thickness  and  lift  distributions  are  taken  to  be  character- 
ized by  the  small  parameters  't  and  oC  , where  ^ / — 0(<  ) . These  parameters 

directly  affect  the  orders  of  the  velocity  components  U*  and  <v ] , in  a 

manner  discussed  below.  The  order  of  the  radial  component,  tT  , is  affected 
only  indirectly  by  the  parameters  and  <X  , in  the  sense  that  the  continuity 
requirements  will  cause  the  radial  velocity  component  to  respond  to  the  per- 
turbations in  the  other  two  components.  In  a properly  designed  blade  row, 
this  overall  continuity  requirement  is  reflected  in  an  axial  rate  of  change  of 
the  annular  cross  section,  whose  magnitude  is  influenced  by  the  thickness  and 
loading  parameters.  For  the  purposes  of  the  present  analysis,  the  rate  of 
change  of  the  annular  cross-sectional  area  is  not  introduced  explicitly  as  a 
small  parameter.  Rather,  the  radial  perturbation  velocity  is  carried  along  in 
the  equations,  with  a view  to  relating  its  order,  through  the  other  two  com- 
ponents, to  the  thickness  and  loading  parameters.  Thus,  the  basic  ordering 
problem  is  treated  in  a quasi  two-dimensional  manner,  involving  only  and 
l/Ly\  # 
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The  appropriate  equations  are  Equation  (28)  and  the  first  component  of 
the  irrotational ity  condition,  here  rewritten  with  the  s - and  vn  -coordinates 
made  dimensionless  by  the  local  chord  length  JC-  J£.\ 

l'-’*'  w.  1 o + *■  sW) 

= (r.,^  ^ ^ 


2 (*V*A>\ 


2 (**/vJ0) 
2(t/c) 


(57) 


The  next  step  is  to  establish  the  dependence  of  the  velocity  components  on  the 
parameter  T . This  dependence  is  influenced  by  the  solidity  of  the  blade  row, 
c / Lt  - At  very  l0W  solidit>r>  the  dependence  approaches  that  of  the  isolated 
airfoil  case.  This  can  be  seen  by  examining  the  blade-surface  boundary  condition, 

a*  = fta , v$'(vc.) 


(58) 


where  A is  of  unit  order.  If  the  equations  above  are  divided  by  and 
A5  * respectively,  where  ~ I ~ ^ follows  that  the  solution  has  the 

familiar  transonic  similitude: 


/ \fJ0  / vJ  a 
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A* 


- 
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(59) 


where 


[_  G J 


(A  = 


I - 


= o (0 


(60) 


and  where  all  the  arguments  appearing  in  the  functions  are  of  unit  order.  The 
dependences  of  and  display  the  well-known  tendency  of  disturbances  in 
two-dimensional  transonic  flows  to  be  felt  at  large  transverse  distances.  The 
full  implications  of  this  similitude  have  been  studied,  both  theoretically  and 
experimentally,  by  Amecke.35 

The  requirement  that  ° be  of  unit  order  cannot  be  satisfied  if 

the  solidity  becomes  of  unit  order,  since  then  A / c is  restricted  to  values  of 

order  one,  while  -t*  O • In  this  range,  a new  phenomenon  comes  into  play, 

namely  that  the  flow  between  the  blades  begins  to  act  like  a transonic  channel 

flow,  in  which  the  distance  over  which  the  transition  through  sonic  conditions 

occurs  is  controlled  by  the  effective  radius  of  curvature  of  the  channel. 

This  influence  was  first  pointed  out  by  Ackeret  and  Rott  . The  implications 

37 

on  the  ordering  of  the  flow  variables  has  been  studied  by  Hall  , who  observed 
that  the  axial  extent  of  the  transonic  zone  (which  affects  the  ordering  of  the 
streamwise  coordinate)  is  governed  by  the  radius  of  curvature  of  the  channel 
throat . 


The  adaptation  of  Hall's  analysis  to  the  present  problem  leads  to  the 
following  conclusions.  The  normal  velocity  component  at  the  blade  surface  is 
still  given  by  Equation  (58),  but  it  is  no  longer  true  that  is  of 

order  ^ , since  the  surface  slope  function  itself  vanishes  at  the  sonic 
location.  In  order  to  establish  the  orders,  it  is  necessary  to  use  the 
irrotationality  condition: 

^ (w-v  VO  I , u 


- r - 


i 


o* ) 


(61) 


where  surface  curvature  is  of  unit  order,  and  does  not  vanish 

in  the  region  of  near-sonic  flow.  Since  v\  ~ O , it  can  be  concluded 
that  a I c " 0 V . In  order  to  deduce  the  orders  of  and  5 c , 

Equation  (57)  is  differentiated  with  respect  to  s , and  the  irrotationality 
condition  is  used  to  eliminate  ■ 

r . >'i,  t.3V  d _ / 

A - _ (Y+'\N\,  + — ' = O t 

I/"*  /C)' 


jr  rs 


(62) 
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From  the  known  orders  of  ^5 / ,jQ  and  , it  is  clear  that  we  must 

have 

Ac  =<=<>)  , *r--C  O’'*) 


The  solution  is  then  of  the  form 


j'J, 


tJ9  , ( _JL- 

i * rx  ' > c.  ) y»  i_i 


£ ' 


(63) 


where  all  of  the  parameters  are  of  unit  order.  It  can  now  be  verified  that 
the  nonlinear  terms  on  the  right-hand  side  of  Equation  (62)  are  smaller  by  a 
factor  ^ than  those  retained. 


It  should  be  noted  that  the  similarity  parameter  ^//Sj  appears  in  a 

number  of  papers  that  deal  with  transonic  channel  flow.  As  an  example,  it  is 

the  appropriate  parameter  in  the  low-frequency  limit  of  Adamson's  study  of 

38 

time-dependent  transonic  channel  flow. 


The  orders  of  magnitude  of  the  radial  velocity  and  the  radial  coordinate 
can  now  be  related  to  those  found  above,  by  using  the  two  components  of  the 
irrotationality  condition  in  which  lF  appears.  The  result,  for  any  value  of 
the  solidity,  is  that  iT  and  F take  on  the  same  orders  as  those  of  tv.,  and 
Vh  . Thus,  the  terms  in  the  small -disturbance  nonlinear  potential  equation 
that  permit  radial  communication  may  be  retained,  in  a consistent  ordering 
scheme . 


It  should  be  stressed  that  the  order-of -magnitude  analysis  above  is 

only  aimed  at  assuring  the  consistency  of  the  potential  equation  derived  earlier, 

which  is  used  here  in  the  sense  of  a composite  equation.  It  is  not  intended 

to  be  a definitive  similitude  study.  Because  of  the  complexity  of  the  problem, 

and  the  large  number  of  parameters  present,  it  is  likely  that  futher  analysis 

will  lead  to  more  precise  definitions  of  similarity  parameters,  and  to  more 

refined  approximations  that  may  be  used  in  the  appropriate  ranges  of  these 

parameters.  Similitude  studies  of  three-dimensional  flows  over  isolated  wings 
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and  bodies  have  yielded  a wide  variety  of  such  parameters.  It  seems  probable 
that  a comparable  variety  remains  to  be  found  for  three-dimensional  turbo- 
machinery flows. 
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CONCLUDING  REMARKS 

The  derivation  of  the  small-disturbance  equations  presented  above  is 
the  first  step  in  a study  of  the  combined  effects  of  three-dimensionality  and 
nonlinearity  in  the  transonic  flow  through  a compressor  blade  row.  The 
high  degree  of  similarity  between  these  equations  and  those  appropriate  to 
isolated  airfoils  in  transonic  flow  suggests  that  many  of  the  numerical 
techniques  developed  for  calculating  the  latter  flows  can  be  carried  over  to 
the  case  of  turbomachinery  flows.  This  carryover  requires  certain 
adaptations,  to  account  for  the  periodicity  conditions  and  the  somewhat 
different  statement  of  the  blade-surface  boundary  conditions.  These  adapta- 
tions have  been  made  succsssfully,  and  are  reported  in  Reference  40. 

The  fact  that  these  adaptations  can  be  made  for  the  small-disturbance 
equations,  suggests  that  it  should  be  a relatively  straightforward  matter  to 
extend  the  adaptations  farther  into  the  nonlinear  range. 

The  equations  derived  above  are  also  useful  in  assessing  the  approxi- 
mations inherent  in  the  use  of  a rectilinear  shear  flow  as  a model  of  a compressor 
flow.  There  are  significant  differences  between  these  models  and  the 
equations  appropriate  to  a blade  row,  especially  in  regard  to  the  representa- 
tion of  three-dimensional  effects.  Fortunately,  the  equations  derived  here 
appear  to  be  no  more  complicated  than  these  qualitative  models. 
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